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Abstract 

It is generally believed that in spatial dimension d > 1 the leading contribution to the entangle- 
ment entropy 5 = —trpA log pA scales as the area of the boundary of subsystem A. The coefficient 
of this "area law" is non-universal. However, in the neighbourhood of a quantum critical point S 
is believed to possess subleading universal corrections. In the present work, we study the entan- 
glement entropy in the quantum 0{N) model in 1 < d < 3. We use an expansion in e = 3 — d to 
evaluate (i) the universal geometric correction to S for an infinite cylinder divided along a circular 
boundary; (ii) the universal correction to S due to a finite correlation length. Both corrections are 
different at the Wilson-Fisher and Gaussian fixed points, and the e ^ limit of the Wilson-Fisher 
fixed point is distinct from the Gaussian fixed point. In addition, we compute the correlation length 
correction to the Renyi entropy 5„ = log trp^ in e and large- expansions. For N ^ oo, this 
correction generally scales as A^^ rather than the naively expected A^. Moreover, the Renyi entropy 
has a phase transition as a function of n for d close to 3. 
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I. INTRODUCTION 



One of the most fascinating and counterintuitive properties of a quantum system is the 
entanglement of its many-body wave-function. In recent years, there has been a lot of interest 
in using entanglement as a theoretical probe of ground state correlations. It is hoped that 
this viewpoint will be particularly fruitful in studying quantum critical points, which realize 
some of the most non-classical, entangled states of matter. 

A useful measure of entanglement is given by the entanglement entropy S, also known as 
von- Neumann entropy. To compute S, we divide the system into two parts, A and B, and 
determine the reduced density matrix pA = trsp, where p is the full density matrix of the 
system. Then, the entanglement entropy, 

Sa = -trApA'^ogpA (1.1) 

If the system is in a pure state, then the entanglement entropy is "mutual", i.e. Sa = SB- 
One may ask how does the entanglement entropy behave near a quantum critical point. 
This question has been addressed completely for one-dimensional critical points with dy- 
namical critical exponent z = 1. Such critical points are described by 1 + 1 dimensional 
conformal field theories (CFT's). In these systems if A is chosen to be a segment of length 
/ and B - its complement in the real line, the entanglement entropy is given by,-i^ 

S=^\ogl/a (1.2) 

where a is the short- distance cut-off and the constant c, known as the central charge, is 
a fundamental property of the CFT. Moreover, if the system is perturbed away from the 
critical point, the entanglement entropy becomes, 

S = A%gaa (1.3) 
o 

where ^ is the correlation length and A is the number of boundary points of the region A. 
Here it is assumed that A and B are composed of intervals whose length is much larger than 

The study of entanglement entropy at quantum critical points in dimension d > 1 has 
received much less attention. The leading contribution to S is believed to satisfy the "area 
law",^ 

S = C-A (1.4) 

where A is the length/area of the boundary between the regions A and B. Physically, the 
area law implies that the entanglement in d > 1 is local to the boundary even at the critical 
point. The coefficient C entering the area law is sensitive to the short distance cut-off, 
and is, therefore, non- universal. So, in contrast to the one-dimensional case, the leading 
term (11. 4p in the entanglement entropy in higher dimensions cannot be used to characterize 
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FIG. 1: The cylindrical geometry considered in calculation of finite size correction to the entangle- 
ment entropy. 



various critical points. 

However, one may subtract the leading non-universal area-law contribution to the entan- 
glement entropy and consider, 

AS = S -C^ (1.5) 

At least for Lorentz-invariant theories that we study here, it is generally believed that 
if d = 2 and the boundary between the regions A and B is closed and smooth, AS* is 
universal. (Additional logarithmic divergences are believed to occur when the boundary 
contains corners/endpoints.^'-) In particular, precisely at the critical point, AS* is just a 
geometric constant. Moreover, AS" is expected to remain universal when the theory is 
perturbed away from the critical point by a finite correlation length ^. 

We note that the above considerations have only been verified by explicit field theoretic 
calculations in free theories. These assertions were also confirmed in strongly coupled su- 
persymmetric gauge theories using the AdS/CFT correspondence .-i^ Recently, Hsu et al- 
found universal corrections for a special class of quantum critical points in d = 2 which are 
described by dimensional reduction to a classical d = 2 field theory. However, such criti- 
cal points are non-generic, and unstable^ in physical situations to quantum critical points 
described by interacting field theories in 3 space-time dimensions. 

In the present work, we compute the geometric and correlation length corrections to 
the entanglement entropy in the simplest generic interacting CFT in d = 2 dimensions - 
the 0{N) model. We verify that these corrections are, indeed, universal. We perform our 
calculations using expansions in e = 3 — c? and Note that the universality of AS* 

formally extends to the range 2 < D < 4, where D = d + 1 is the space-time dimension.^ 

In the rest of this paper we consider the following geometry. We take two semi-infinite 
regions A and B with a straight boundary at x = 0. The boundary extends along the 
remaining d — 1 spatial directions, each taken to have a length L. For technical reasons, we 
impose anti-periodic boundary conditions along each of these directions. We also consider 
more general boundary conditions with a twist by an arbitrary phase in a theory of N/2 
complex scalar fields. So in the physical case d = 2, our space is an infinite cylinder divided 
into regions A and B along a circle of length L, see Fig. [H In this geometry the entanglement 



In D = 4, S develops new singularities associated with the extrinsic curvature of the boundaryii 
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entropy at the critical point is given by, 



S = C^ + 1 (1.6) 

We explicitly compute the universal geometric constant 7. To leading order in 4— e expansion 
we obtain, 



Ne 

7 -- 



(log 6*1(^^—^—^,2) — y log 77(2) ) , D = 4— e, Wilson-Fisher fixed point 

\ 2n 47r / 

(1.7) 

Here 61 and r] are Jacobi elliptic and Dedekind-eta functions. The sign of 7 depends on the 
value of ip: it is negative for = vr (anti-periodic boundary conditions) and positive for 
9? — > 0. Note that eq. (11.71) is only valid for ip ^ e^/^. For zero twist (periodic boundary 
conditions), we hypothesize that to leading order, 

Ne 

The result (11. 7p should be compared to the corresponding value at the Gaussian fixed point 
in 4 — e dimensions. 



(log 6i(— — —— log rj (i)] , D = 4 — e, Gaussian fixed point (1.9) 

\ 2n An J 



We see that I7I is parametrically smaller at the Wilson-Fisher fixed point than at the Gaus- 
sian fixed point. Thus, entanglement entropy distinguishes these two fixed points already 
at leading order in e expansion. 

If we perturb the system away from the critical point, we can take the limit L ^ 00 and 
obtain the general scaling relation, 

Td-l jd-l 

S = C^ + r^ (1-10) 

where r is a universal coefficient that we compute. In general, one has to make a specific 
choice for the definition of the correlation length ^. In the 0{N) model there is a very 
natural choice, ^ = m~^, where m is the gap to the first excitation. Note that in the present 
work we only consider the phase of the 0{N) model with unbroken symmetry. The value of 
r to leading order in 4 — e expansion is found to be, 

N 

D = 4 — e, Wilson-Fisher fixed point (l-H) 



1447r' 

As with the finite size correction, r is parametrically smaller at the Wilson-Fisher fixed point 
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than at the Gaussian fixed point where,- 



, D = 4 — e, Gaussian fixed point (1-12) 



247re 

In addition to the entanglement entropy, we study the Renyi entropy, 

Sn = -^\ogtrAP^A (1-13) 
1 — n 

The Renyi entropy always naturally appears in field-theoretic calculations as it is related 
to the partition function of the theory on an n-sheeted Riemann surface. One then obtains 
the entanglement entropy by taking the limit, S = lim^^i Sn- At least for n close to 1, 
the Renyi entropy is believed to possess the same universal properties as the entanglement 
entropy. In particular, the finite size and correlation length corrections are given by, 

Sn = Cn^+ln (1.14) 
Td-l Td-l 

Sn = Cn^ + r^^ (1.15) 

where the non-universal coefficient C„ of the leading area law term, as well as the universal 
coefficients 7„, r„ are now n dependent. We compute r„ in 4 — e and large- expansions. 
A careful renormalization group analysis demonstrates that r„ is parametrically enhanced 
in both of these limits. In particular, r„ ~ 0{\) in the 4 — e expansion. However, the 
enhancement is most striking in the large- expansion where we find ~ 0{N'^). Such 
scaling is in contrast with the result r„ ~ 0{N) that one would obtain at each order in 
for fixed correlation length ^, implying that the limits ^ — >■ cx) and — >^ cx) do not commute. 
As far as we know, this is the first violation of naive large- A^ counting in the 0{N) model. 
A common feature of the two expansions is that the leading term of r„ behaves as r„ ~ n — 1 
for n 1 and does not contribute to the entanglement entropy S. Hence, r ~ 0{N) in the 
large A^ limit and r ~ 0(1) in the 4 — e expansion. 

Another unusual phenomenon that we find in 4 — e expansion is non-analytic dependence 
of the coefficients 7„, r„ on n. In fact, 7„ and r„ will have a discontinuity at n = n*, where n* 
is generally non-universal and lies in the range, 1 < n* < l + |^±|e. The n-dependence of 7„ 
and r„ for n < n* and n> n* is, however, universal. Thus, we have two universal branches 
for 7„ and r„. We note that eqs. f ll.l4p and f ll.lSp are understood in the limit L — > cx), 
—>■ oo. However, there appears a new divergent length-scale in the problem as n — > n*, 
and the limits n ^ n* and L — ^ oo, ^ — > oo do not commute. In particular, if we fix the size 
of our regions L or the correlation length ^, the n-dependence of the Renyi entropy Sn will 
be completely analytic. Moreover, due to the emergence of a new length-scale as n — > n*, 
in the crossover region Sn is not entirely universal. We stress that any non-analyticity and 
non- universality only occurs away from the point n = 1. In particular, the entanglement 
entropy S = lim^.,! Sn is well defined and universal. 
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The non-analytic behaviour discussed above is also found to occur in the large- expan- 
sion in dimensions 2.74 < d < 3. The limited range of d suggests that this phenomenon 
might be absent in the 0{N) model in the physically relevant case d = 2. Nevertheless, 
we expect that such non-trivial n dependence will occur quite generically at other quantum 
critical points. 

This paper is organized as follows. In section [Tll we remind the reader of the replica trick, 
which relates the entanglement entropy to the partition function on an n-sheeted Riemann 
surface. In section llllt we show that the coefficient of the correlation length correction to 
the Renyi entropy r„ is parametrically enhanced in both expansions we consider. Sections 
IIVI and |V] are respectively devoted to the evaluation of correlation length and finite size 
corrections in 4 — e expansion. In section |Vl] we compute the coefficient r„ in the large- 
expansion. Some concluding remarks are given in section IVIIi 



II. THE REPLICA TRICK 

We consider the 0{N) model in D = d + 1 space-time dimensions. The action for the 
iV-component real scalar field is given by , 



S = I d^'xdr ( -{d^4>f + + ^4,^ ] (2.1) 



We divide our space into two regions A and B with the boundary being a d — 1 dimensional 
plane at x = 0. We will denote the coordinates along the boundary directions by x^. The 
Renyi entropy Sn may be calculated as, 

Sn = j^\og^ (2.2) 
1 — n z " 

from which we obtain the entanglement entropy, 

S = lim Sn (2.3) 

Here is the partition function of the theory on an n-sheeted Riemann surface. This 
Riemann surface lies in the a;|| = (t, x) plane and has a conical singularity at (r, x) = (0,0). 
The surface is invariant under translations along the x± directions. We may use the following 
metric for our space-time, 

ds"^ = dr'^ + r'^de'^ + dx\ (2.4) 

where r, 9 are the polar coordinates in the (r, x) plane. Concentrating on this plane, we see 
that the metric is exactly the same as for the usual Euclidean plane; the only modification 
is that the angular variable 9 has a period 9 9 + 27m. 
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III. PARAMETRIC ENHANCEMENT OF CORRELATION LENGTH CORREC- 
TION 



In this section, we show that the coefficient r„ of the correlation length correction to the 
Renyi entropy, eq. (II. 151) . is parametrically enhanced in both expansions that we consider. 
Moreover, we demonstrate that r„ can to leading order be extracted from the properties of 
the theory at the critical point. 

We start with the 0{N) model perturbed away from the critical point t = tc hj a finite 
i = t — tc > (we drop the tilde below). To compute r„, we need to find the dependence of 
the partition function Z„ on the mass gap m = Here we assume that the dimensions 
of the boundary L ^ ^, so that we can take the limit L ^ oo. It is useful to differentiate, 

tM% = d''x{<p\x))^ -n I rf^x(0^(x))i) (3.1) 

■^l ^ V^n-sheets ^1-sheet / 

^L'-^ I d^x.{{<p\x))n-{<P\x)),) (3.2) 
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n— sheets 



where we have used the fact that the contribution to the integral from each of the sheets is 
the same (from here on, all integrals over d^x\\ are understood to be over ra-sheets). Now, 
recalling, m f^, we may convert the derivative with respect to t into a derivative with 
respect to m, 

""i^^^ J = d'xiim'ix))n-{^'ix)),) (3.3) 

The expression t((0^(x))„ — (0^(x))i) is renormalization group invariant.^ Thus, we may 
write, 

t((02(x))„ - (</)^(a;))i) = m^/„(mr) (3.4) 

where /„ is a universal function. The function /„ is expected to decay exponentially for 
mr 3> 1, and the integral in (13. 3p converges for r — > oo. The short- distance asymptotic 
of fn is controlled by the critical point. From the scaling dimension of the operator 0^(x), 
[(^^(x)] = D — u'^, we conclude, 

where dn is a universal constant. So the integral in (13.31) converges for r — * 0, provided that 
z/^^ > D — 2.^ In the 0{N) model in both expansions we consider, z/^^ = D — 2 + z/i, where 



^ Two subtractions (constant and linear in i), in addition to the multiplicative renormalization, are needed 
to render the operator 0^ finite. However, these subtractions cancel among the two expectation values in 

^ Otherwise, a UV divergence appears which adds a piece analytic in t to the entanglement entropy, in 
addition to the singular contributions discussed below. 



7 



the correction ui is given to leading order by, 



ui = D = 4-e (3.6) 

1 ST{D) / _ X _ 32 

- iVDr(2-D/2)r(D/2-i)2r(D/2)' ^^^^"^^"s;^' ^^'^^ 

In particular, i/i > and i^"^ asymptotically approaches D — 2 from above in both limits. 
With these remarks in mind, we integrate eq. (13.31) with respect to m, 

log^(t)-log^(t = 0) = -^^^^(mL)^-^y^ rf« «/„(«) (3.8) 

This is as far as we can proceed in general - to make further progress one needs the function 
fn{u). However, we have already noted that due to the fact, u^^ ^ D — 2, the integral in 
(13. 8p is very close to diverging in both expansions. Hence, to leading order in e or this 
integral is saturated at short distances. 
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and 

log§^-— rf„(mL)'^-i (3.10) 

where we've dropped the constant contribution at the critical point t = 0. So, the universal 
coefficient r„ of the correlation length correction, eq. (ll.lSp . is given by, 

'''' -rfn (3.11) 



(1 - njui 

Thus, to leading order the problem is reduced to evaluating the coefficient dn in (13.51) . Since 
this coefficient is a short distance property, we may work directly at the critical point. Note 
in particular that in the large N limit, dn ~ 0{N), so our result for log scales as iV^. 
This is in contrast to the linear in N behaviour that one would obtain at any finite order in 
the 1/A^ expansion for a fixed correlation length ^. 

It turns out that the leading term (I3.1ip behaves as r„ ~ [n — 1) for n ^ 1 in both 
expansions and does not contribute to the entanglement entropy, eq. (12. 3p . Thus, the 
correlation length correction to the entanglement entropy has the expected scaling r ~ 0{N). 
To proceed systematically beyond the leading order one needs to use renormalization group 
(RG) technology that will be developed explicitly in the context of 4 — e expansion in section 

iivcni 
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IV. 4 - e EXPANSION: CORRELATION LENGTH CORRECTION 



In this section we compute the correlation length correction to the entanglement entropy 
in 4 — e expansion. Recall that for the interacting 0{N) model, ui = — (D — 2) ~ 0{e) 
m D = A — e dimensions, hence the argument in section IIIII can be applied. This is also true 
for the non-interacting (Gaussian) fixed point for D = 4 — e, where ui = e, allowing us to 
compare the predictions of our method to the exact calculations of Ref. 0. We first consider 
the Gaussian fixed point and then proceed to the Wilson-Fisher fixed point. 

A. Gaussian theory 

Consider the Gaussian theory, 

L = l{d,<l>r + 1<I>' (4.1) 

where, t = m^. We need to compute the expectation value, 

(02(x)). - (02(x))i (4.2) 

at the critical point, t = 0. To leading order we may work mD = A. The massless propagator 
on an n-sheeted Riemann-surface in D = 4 is known to be,— 

^ / , ^ N sinh(r7/n) 

Gn r, r', 9, x^) = ^ \^ 4.3 

oTT^nrr' smh ?7(cosh(?7/n) — cos[U/n)) 



where 



Hence, 



coshr/ = (4.4) 



(0m-(0^WK = ^(^-l) (4.5) 



So comparing to eqs. (13.41) . (13.51) . we obtain, 

dn = ^7 f ^ "~ 1 I 5 Gaussian fixed point, D = A — e (4.6) 

487r"' J 

We can now use eq. (13.111) to compute the coefficient r„ of the correlation length correc- 
tion. As noted above for the Gaussian theory, vi = e, so 



487re \ n 
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and for the entanglement entropy proper, 



r = lim r„ = (4.8) 

n^i " 247re ^ ^ 

This can be compared to the exact result of Ref. [2I, 

Eq. (14.91) is in agreement with our result (14.71) to leading order in e, which is all that the 
discussion in section |TTl] guarantees. 



B. Interacting theory 

We now proceed to consider the interacting 0{N) model, eq. ( 12. ip . We again need to 
compute the expectation value (14. 2p . Naively, one would expect that at leading order in e, 
one can work with the mean-field approximation, u = 0, recovering the result (14.61) . Then, 
one would simply substitute (14.61) into eq. (13. lip and use the appropriate z/i, eq. (13. 6p . for 
the Wilson-Fisher fixed point. However, such reasoning turns out to be too simple minded, 
as it neglects "boundary perturbations." Indeed, our conical singularity will generally induce 
local perturbations at r = 0. Of these, the term with the lowest engineering dimension is, 

SS=^j d^-^x± 02(r = 0, x±) (4.10) 

In the absence of the conical singularity this perturbation is known to be irrelevant in the 
0{N) model as the scaling dimension [c] = z/^^ — 2 < 0.— However, as we will now show, 
the presence of the conical singularity will modify the renormalization group flow of the 
coefficient c. 




FIG. 2: Leading correction to the propagator 5Q ' due to the boundary perturbation. Here and 
below, a cross denotes an interaction vertex of c. 



The engineering dimension of the coupling constant c is zero in any space-time dimension 
D. We wish to compute the /3-function, /3(c). Let us perform perturbation theory in u and 
c for the two-point function {(j)a{x)(l)p{x')) = 6a/3G{x,x'). It is sufficient to work in D = 4 
dimensions to compute the leading terms in /3(c). We use a mixed momentum/position p^, 
x\\ representation. To first order in c and zeroth order in u, we have the simple diagram in 
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y X 

a) 



R S E 
X » X ^ 



b) 



FIG. 3: Corrections to the propagator, a) 6Q^'^ and b) 6Q'^'^. Here and below, a dot denotes an 
interaction vertex of u. 




FIG. 4: Corrections to the propagator 6Q 



Fig. El 

= -cG„(x|,,0,px)G„(0,xj|,p^) (4.11) 

where the superscripts on 6 indicate the order in c and u. Notice that the bare propagator 
Gn{x, x'), eq. (14.31) . remains finite as its arguments approach the conical singularity. In fact, 

GniO,x) = -G,ix) (4.12) 
n 

Also, Gn{x\\,x'pp±) is just the two dimensional massive propagator (— V| + Pj_)~^ on an 
n-sheeted Riemann surface. In particular, Gn{xi\,0,p±) = ^Kq{p±\x\i\) (which implies that 
the relation (14.121) is actually correct in any dimension). Thus, the correction (14.111) is finite. 
We next consider the Hartree-Fock (first order in u) correction to the propagator. Fig. [3] 

a), 

5°''^g{x\\,x'pp±) = -{N + 2)u j d^y\\ Gnixi\,y\\,p±)Gniy\\, x\^,p±){Gniy,y) - Gi{y,y)) 

(4.13) 

We have already evaluated Gn{y, y) — Gi{y, y) ~ -y, eq. (14. 5p . Thus, the integral (14.131) has 



11 



an ultraviolet divergence in the region y\\ 0, 

5'Mh,x[.Pa-) = - ^) G'„(x|i,0,p^)G„(0,4,p^)log(A) (4.14) 

Notice that this divergence is local to the conical singularity and, as is evident from eq. 
(14. lip , can be canceled by an additive renormalization of the coupling constant c. Hence, the 
perturbation fl4.10p will be automatically induced by the presence of the conical singularity. 
We also consider the second order contribution in c to the propagator. Fig. [3]b), 

52.0^(0:11, 4, px) = c2G'„(a;||,0,p^)G„(0,xj|,p^)G„(0,0,p^) (4.15) 

The quantity G'n(0,0,p_L) is UV singular. 



G„(0,0,px) = j rfVG'„(0,0,i/^)e-^^^^^^ = d'y±^e^P^y^ iog(A/pj 

so 

^2 



(4.16) 



52.0^(0:11, 4, j9^) -^G„(x|,,0,px)G„(0,xj|,px)log(A) (4.17) 

II 27m " 



The divergence of fl4.17p is a manifestation of the well-known fact that the two-dimensional S- 
function potential requires regularization. Again, from (14. lip , we observe that the divergence 
can be eliminated by a renormalization of the coefficient c. 

Finally, we consider corrections which are bilinear in c and u, Fig. HI For c - small, these 
corrections are generally subleading compared to Fig. [3] a). However, for n — >■ 1, 

vanishes, and the diagram in Fig. H] c) becomes important. On the other hand, the 
diagrams in Figs. H] a,b) can be ignored to leading order for all n since they also vanish at 
n = 1.^ With this in mind, we only need to evaluate Fig. |l]c) at n = 1. We recognize, that 
this is just the diagram corresponding to the usual multiplicative renormalization of the 0^ 
operator. Explicitly, 

5^'^g{x\i,x'pp±) {N + 2)uc j £yi\Gi{x\i,yii,p±)Gi{y\\,x'pp±) j d'^z±Gi{y\i, z±f 

1 



= {N + 2)uc j d y\\Gi{x\\,ypp^_)Gi{yi\,x'pp±j ^^^^^^ 
UV {N + 2)uc 



87r2 



y\\ 

G'i(x||,0,px)Gi(0,x;|,px)logA (4.18) 



Technically, these diagrams contain (log A)^ divergences, and one needs to use a consistent regularization 
method to evaluate them. 
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We can now introduce counterterms to cancel the divergences considered above, 

where Cr and are the renormahzed couphng constants and /i is the renormahzation scale. 
Note that the coefficient of the u^Cr term has been only computed at n = 1. So, 

/J(c,) = ,|-c.[__ = '^:.!>"- f „ - + <^!^±|^ + ± (4.20) 



247r \ n Stt^ 27m 



Note that the RG flow of u is not affected by the boundary perturbation or by the presence 
of the conical singularity, 

I3{ur) = -eur + ^-^^ui (4.21) 



■in 

and we have the usual Wilson-Fisher fixed point u* 



iV + 8 2 

r 





FIG. 5: /3-function of the boundary coupling for a) Non-interacting theory (u = 0), b) Interacting 
theory, n = 1, c) Interacting theory, n < ric, d) Interacting theory, n > ric- 



We now discuss the RG flow of Cr in detail. Let us start with the non-interacting theory, 
u = 0, which corresponds to the well-studied problem of a particle in a two-dimensional 
5-function potential. Then, P{cr) = 2^c^- As demonstrated in Fig. [5] a), the coupling 
constant flows logarithmically to zero for > and runs away to — oo for Cr < 0, 
signaling the formation of a bound state. 
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Next, consider turning on the interaction u, in the absence of conical singularity {n = 1). 

2 

Then, l3{cr) = —r]2{ur)cr + where 7]2 is just the usual anomalous dimension of the 0^ 
operator, ([0^] = D — 2 — ri2), 

, , (N + 2)Ur 

V2iUr) = ^ , (4.22) 

The RG flow of c is sketched in Fig. [5]b). We find two fixed-points: c,^ = and c~ = 
— ^i|(27re). The first fixed point c+ = is stable, due to P'{cr = 0) = —rj2{u*) > 0, which 
implies that for c - small, the perturbation fl4.10p is irrelevant.^- This conclusion can be 
immediately reached by consideration of scaling dimensions at the interacting fixed point, 
since [c] = D - 2 - [0^] = < 0. 

The second fixed point c~ is unstable, and for < c~ the RG flow runs away to = —00. 
Naively, such a flow may be interpreted as a tendency of to condense in the vicinity of r = 0. 
However, this would result in a condensate that is effectively D — 2 <2 dimensional, which, 
at least for N > 2 and t > 0, is prohibited by the Mermin- Wagner theorem. Exactly at 
the critical point, long-range forces could, in principle, stabilize the condensate. However, 
as we will discuss in section IVIt large- iV expansion suggests that no such condensation 
occurs even at t = 0, and the flow actually terminates at a scale invariant fixed-point, 
which is inaccessible in our perturbative expansion. However, this fixed point can likely be 
interpreted in terms of a fluctuating "boundary" order parameter. 

Finally, we proceed to the interacting case in the presence of a conical singularity. For 
n < nc ~ 1 + f we again obtain two fixed points. Fig. [5]c), 



= " I -NTS'"' ^ V [nTS ) " ' - 3 ivTi^" ' ^'-''^ 

The fixed point is stable, while c~ is unstable. In the limit n — > 1, which is relevant for 
the computation of entanglement entropy, smoothly evolves to the = stable fixed 
point, which we obtained in the absence of the conical singularity. Moreover, for n — 1, we 
expect the starting point of the RG flow c,. 0. Hence, for n close to 1 the RG flow will 
terminate at the fixed point c^. Thus, the main effect of the conical singularity is to shift 
c,^ away from 0. The parametric magnitude of this shift depends on whether 1 — n ^ e or 
11 - nl < e: 



2 N + 2 

, 27r, , 27riV + 8(ra-l)2 , , ^ ^ 

"-T<"-"-y 7V + 2 e • <^-2^' 

Thus, for 1 — n ^ e, ~ 0{^/e)^. this is the regime in which the UrCr term in the /3-function 

fl4.20p can be ignored. On the other hand, for |n — 1| <^ e, ~ (n — 1) <^ e and the UrCr 

term in l3{cr) becomes important. Note that in both regimes, is parametrically small and 
the perturbative expansion in Cr is justified. 
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For n > He, both fixed points disappear and the RG fiow runs away to c,. = — oo, Fig. [5] 
d). As discussed above for the case n = 1, large analysis suggest that the fiow is towards 
another fixed point (which itself evolves as a function of n). Now there are two possibilities. 
If as n increases from 1 to ric, the initial value of c^, determined by the microscopic details 
of the theory, satisfies Cr{n) > c~{n) then the run-off to the Cr = — oo fixed point will occur 
precisely at n = n* = Uc- On the other hand, if the initial value of the coupling Cr{n) < c~{n) 
for n > n* where 1 < n* < ric, the runaway to Cr = — oo will occur before n reaches tIc- Note 
that the value of n* is generally non- universal. In either case, the long-distance physics is 
controlled by the fixed point for n < n* and the Cr = — oo fixed point for n > n*. Thus, 
the constants 7„, r„, eqs. (11.141) . (ll.lSp will always have a discontinuity at some n = n*, 
1 < n* < Tic- Note that eqs. (]1.14p . ( Il.isp are understood in the limit when the size of the 
regions whose entanglement entropy we are computing and the correlation length tend to 
infinity. However, as n — > n* a new divergent length scale emerges in the problem. In fact, 
we can think of the point n = n*, t = as a multicritical point. Thus, the limits L, ^ — > oo 
and n ^ n* do not commute. In particular, if we fix L or ^, the dependence of the Renyi 
entropy on n will be completely analytic. Moreover, the emergence of a new length-scale as 
n n* implies that the Renyi entropy in the cross-over region is not entirely universal. 



FIG. 6: Leading contributions to {(jP'{x))n (denoted by a black square here and below): a) Mean- 
field result, b) Correction due to the boundary perturbation. 

Having discussed the non-trivial n-dependence of the Renyi entropy that occurs for n 
away from 1, we come back to the range n < n^. and concentrate on the fixed point. We 
will from here on denote as c*. Let us now compute the value of {(fP'{x)) at this fixed 
point. The leading correction to the mean-field result. Fig. [6] a), eq. (14. Sp . is given by the 
diagram in Fig. [6]b), 




a) 



b) 




(4.26) 



Since to leading order we still have t = m 



2, from eqs. ([32D and fl33D . 




(4.27) 



15 



and from eqs. (13. 6p . (13.110 . the coefficient of the correlation length correction to the Renyi 
entropy is, 

7m(A^ + 8) , 

rn ^ --r\, fd^ 4.28 

oe(l — n) 

As we see, in the regime 1 — n ^ e, taking the boundary perturbation into account only 
weakly modifies the mean-field result for eq. (14.61) . by a term of order ^/e. Note that r„ 
is still strongly modified due to a different value of Vi. 
However, in the regime |1 — "n-l <^ e, 

N{N + d>){n-lf . 
(iV + 2) 72.^e ' ^'-''^ 

N(N + 8? n - 1 , , , , 

iV + 2 432..^ - 

Thus, for n — s> 1, the behavior of dn at the Wilson-Fisher is drastically different from the 
mean-field result, eq. (14. 6p . In particular, notice that to the present order in e, the correction 
due to the boundary perturbation precisely cancels the term linear in — 1 coming from 
eq. (14.51) . The technical reason for this remarkable cancellation is as follows. For n —>■ 1, 
we expect Cj. ~ 0(n — 1), and we can work just to first order in c. Then, in considering 
the corrections to the propagator, we can drop the diagram in Fig. [3]b), keeping only Figs. 
[3] a) andlUc). These diagrams are, essentially, Hartree-Fock corrections to the propagator, 
and the "Hartree-Fock potential" at y is just {(f)'^{y))n — {4>'^{y))i ~ V^/y- a result, the 
diagrams diverge for 0. The /3-function for the coupling constant Cr vanishes precisely 
when this divergence is absent, i.e. {(t?{y))n — = 0. 

The crucial consequence of eq. (14.291) is that to this order the correction to entanglement 
entropy proper, r = lim„^i r„ = 0. Thus, 

r~0(l), D = A-e (4.31) 

We conclude that the correlation length dependent contribution to the entanglement entropy 
at the Wilson-Fisher fixed point is parametrically smaller than at the Gaussian fixed point 
in D = 4 — e, eq. (14.81) . As a result, we have to proceed to higher order in e to evaluate it. 
This will be done in the next section. 

Before we perform the higher order computation, let us ask how do the correlation func- 
tions of the field 0(x) behave as x approaches the conical singularity. This question is 
connected to the effective boundary conditions on the field cf) that are generated at the sin- 
gularity. In accordance with the general theory of boundary critical phenomena,— we expect 
the field to satisfy the operator product expansion (OPE), 

(t){x\\,XA_) ~ r"(/)(0, xx), r ^ (4.32) 

where 0(0, x^) is an operator living on the conical singularity. The exponent a can be 
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extracted from the two-point function Q{x,x'). Combining the free propagator with the 
boundary correction, eq. (14.111) . 



g{xii , 4 , p±) = ° f 1 + log(p_Lr)) G„(0, 4 , p±) (4.33) 
from which we conclude, 

a=^ (4.34) 

27rn ^ ^ 

Note that from eq. fl4.23p the exponent a is positive for n < 1, implying effective Dirichlet 

boundary conditions on (f){x) at the conical singularity. On the other hand, a is negative for 

1 < n < ric and correlation functions of (f){x) exhibit a power-law divergence as approaches 

the origin. 



C. Beyond the leading order in e 

1. The inhomogeneous renormalization group equation 

At leading order in e, our calculation has relied on the integral in eq. (13. 8p being saturated 
at short distances, u = mr — ^ 0, allowing us to work directly at the critical point. However, 
we saw that the coefficient dn of the short-distance asymptotic of /„, eq. (13. 5p . behaved as 
dn ^ {n — l)^/e for n 1, giving no contribution to the entanglement entropy. We expect 
that to next order in e, dn will acquire a term linear in n — 1, (i„ ~ e{n — 1), which by 
eq. (13.1 op will give a contribution of 0(1) to 5*. Notice that this is of the same order as 
the contribution of the long distance, u oo, part of the integral (13. Sp . which now has to 
be taken into account. Thus, we need to compute the long distance part of /„ to leading 
order in e and the short distance part to subleading order. Although the separation between 
short and long distance contributions is unambiguous to present order, it is convenient to 
introduce a formalism that allows one to consistently treat the problem order by order in e.^ 

Let us define, 

m =n [ rf^xii (( mx)]^)^ - { mx)]^)^) e-^^^' (4.35) 

J 1— sheet 

Here, we have introduced the usual renormalization of the 0^ operator, 

[^\x)]r = ^<P\x), (f ) 't (4.36) 

We are considering $ at a finite momentum p in order to make $ well-defined even at the 



^ We note that the discussion below closely parallels the renormalization group technology used to calculate 
the specific heat in the classical 0{N) model. 
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critical point, t = 0. We are actually interested in computing $ at p = in the gapped 
phase, t 7^ 0, as from eq. fl3.2p . 

t.|;log§ = -^t.$(p = 0)L^-^ (4.37) 

As already observed in section IIIH although the integrand in fl4.35p is finite, the integral 

diverges logarithmically for |x| — > at each order in u. Thus, requires an additive 
renormalization, 

^(p) = ^rip) + CiUr, Cr, yU/A)/i"' (4.38) 

where C is a renormalization constant. We will use dimensional regularization below, so 
that C is, in fact, just a function of Ur and c^. Then satisfies the inhomogeneous renor- 
malization group equation, 

d d d /' (9 \ \ 

with 

f d d \ 

B{ur,Cr) = - [f^Mg;^ + f3{cr)-^ - {^2^ + e) j c,) (4.40) 



where as usual, 



V2M = 



log^ (4.41) 



Note that B must be finite, as the left hand side of eq. f l4.39p is finite. The solution to 
fl4.39p can be represented as a sum of the solution to the homogeneous RG equation and a 
particular solution. In the scaling limit, tr ^ 0, 

MP = 0) = Asfi" i-^] + AnsiUr, Cr)fl-' (4.42) 

where the coefficient of the particular solution A^s satisfies, 

d d \ 

Hence, at the critical point, 

AnsiK. <) = —B. = (4.44) 

r]2 + e ui 

where we recall our definition in section llllt vi = v' -{D - 2) and z/"^ = 2 + t]2. 
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Thus, from eq. fOTIl . 



7, A ( (f \''\^-'^ 



where we've dropped terms analytic in tr- Note that the mass gap m is related to /i 

via a finite proportionality constant, which at leading order in e is just 1. So to leading 

order, 

(4.46) 



2(1 -n) 

Hence, we must compute As. To do so, we perturbatively calculate $r(p = 0) and B{ur, Cr). 
As can then be determined by matching the perturbative expansion with the solution to the 
RG equation 04.421) a the critical point, where the corrections to scaling vanish. Notice that 
we always need to compute B to one higher order in e than $r.(p = 0) due to the factor ui in 
the denominator of eq. (14.441) . Moreover, since $r is finite for e — > 0, while Ans = — -B^,/z/i 
behaves as 1/e, to leading order As = —Ans = B^/ui. Precisely this fact was utilized in 
section [IIT| and we identify to leading order B^ = 2Tmdn- 



2. Regularization 

For the purpose of computing the entanglement entropy S we can work to linear order in 
n — 1. Since the fixed point value c^, ~ 0{n — l), we also work to linear order in c. Therefore, 
all diagrams that include an insertion of c can be evaluated at = 1. In addition, power 
counting indicates that if we work to linear order in c, all diagrams will be finite for D < 4 
(by contrast, higher order diagrams in c, such as Fig. [3]b) diverge even for D < 4). Thus, 
we use dimensional regularization and minimal subtraction below. We remind the reader 
that in dimensional regularization the bare coupling constant u = fi'^Ur Zu/Z^. We hst below 
the renormalization constants in the MS scheme to the order that they will be needed in 
our calculation. 

I, (AT + S) ^ 

Z2 ^ {N + 2)Ur {N + 2){N + 5) / Ur \^ 5{N + 2)/Ur\^ 



Correspondingly, 



= -(AT + 2)^ (4,50) 
As we saw, the boundary coupling constant c will also require renormalization. To linear 
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order in c, 

C = D{Ur) + yC, (4.51) 

where we observe that the muhiphcative renormahzation of c to zeroth order in (^ — 1) is just 
Z2/Z. On the other hand, the additive renormahzation, which behaves as D[ur) ~ {n — 1) 
for n — > 1, needs to be computed exphcitly. So the /3-function, 

/?(c.) = -(y) f^M^ - ^2{Ur)Cr (4.52) 



3. Entanglement entropy to 0(1) 

To calculate the entanglement entropy to 0(1) in e, we need to find the finite part of 
$(]9 = 0), eq. (14.351) . at t 7^ to 0(1) in m and the divergent part of $(p), which determines 
B, eq. flCTj) . to 0{u). 

$(]?) to 0(1) in u is given by the two diagrams in Fig. [61 The diagram Fig. [6] a) is just 
the mean field contribution computed in Ref. 0, 



$(]9 = 0)mf = N J d x\\ {Gn{x, x) — Gi{x, x)) 

= nJ j d\ {G^=\x, X- kl + m')-n^ l) 

N f 1\ f d^-'^ki 1 
n 



12 V nJ J (27r)^-2 kl + rn^ 
iV / l \ r(2-Z}/2) . 4 



12 V (47r)^/2-i ^ ^ 

where G^^^(x, x'; M^) is the two dimensional massive propagator on the n-sheeted Riemann 
surface, and we have used the result proved in Ref. 0, 

j d'x\\ {G^='{x, X- M') - G?=\x, X- M')) = -1 (^n - ^ (4.54) 

The diagram in Fig. [6]b) is the boundary correction, 

5i'0$(p = 0) = -iVc. 1 d^x\\ j d''^%^Gl{x\\,y^) = -Nj-^^-^m''-' (4.55) 

Combining eqs. (14.531) . (14.551) . 

$(p = 0) ^4'^ -N + 1^^) (7 + - i - I^SW/^)) (4-56) 
where we keep only terms linear in n — 1. 
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Subtracting the pole, we obtain for the additive renormahzation constant C, eq. fl4.38p . 



c''^'-n(^ + ^]- (4.57) 



and consequently from eq. (14.401) . 



and 



b''^\c = -n(^ + ^] (4.58) 
V 127r 8n^ ' ^ ^ 



$.(p = 0) -N + ^] Qlog47r -1- log(m//i) ) /i^ (4.59) 



In particular, at the critical point, by eq. (I4.25p . 

* o(i) 27r 

and 

$:(p = 0) = 0(e), 5, = 0(e) (4.61) 

Thus, in the minimal subtraction scheme $*(p = 0) vanishes at the critical point to 0(1) 
in e. The fact that B^, = 2'Kndn vanishes to 0(1) in e has already been observed in section 
HVBl Thus, from eqs. ( 021) . (gaU), 

A. ^ (4.62) 

We now proceed to evaluate B to 0(e). To do this, we compute $(p) at the critical 
point. We first evaluate ([0^]r)n — ([0^]r)i and use it to determine the renormahzation of 
the coupling c in dimensional regularization. We then perform the Fourier transform, eq. 
f l4.35p . to find the subtraction constant C and hence B. To leading order, we have the two 
familiar diagrams in Fig. [6l 



(02(x))„-(02(x))/i='^iV 
where we've defined, 



J{D) - Cr 



T{D/2 - 1 



\3 



167r^/2+ir(Z} - 2) 



y,D-2 



(4.63) 



G„(x,x)-Gi(x,x) = ^ (4.64) 

Note that in dimensional regularization = NGi{x^ a;) = at the critical point. We will 
show in section HV C 41 that to linear order in n — 1, 
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In particular, J{D = 4) = — |^ in agreement with eq. (14.50 . We note that the diagrams 
that contain the tadpole (14.641) can effectively be evaluated with n = 1. The computation 
is simplest in position space, where one uses. 



Gi(x, x') 



T{D/2-l) 



47r^/2| 



t\D-2 



(4.66) 






FIG. 7: Contributions to {4P'{x))n — {4P'{x))i at order u. The counterterm 5^c is denoted by a 
circled cross liere and below. 



At order m, ([0^])n— receives contributions from the diagrams in Fig. [71 Note that 

the diagram c) is the renormalization of the coupling constant cq = Cr + + .... Taking 
the multiplicative renormalization of the operator cf? into account, we obtain. 



^)n-mr)l = N 



Z2 1 V{D/2-l)T{2- D/2f {N + 2)uriJi' 
~ 167r^/2(Z) - 3)F(4 - D) r2(^"3) 



X ( J rp/2-2)^ \ _ F(Z^/2-l)3 ^ 

1 l67r^/2+iF(Z} - 2) 7 167r^/2+iF(D - 2) ^ ^ 



Performing minimal subtraction, 

(^V + 2K/^^A (4.68) 

Notice that the coefficient of the multiplicative renormalization is precisely Z2/Z a,s expected. 
We also obtain the additive renormalization constant, eq. (I4.5ip . 

[N + 2)ur n-1 , , 



Hence, from eq. (14.521) . to first order in u, 



^icr)''^^\N + 2)uJ^ + ^^] (4.70) 
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in agreement with the expression (14.201) obtained earher using cut-off regularization. 

By Fourier transforming eq. (I4.67p . we can compute $(p) at the critical point to order 
u. From the divergent part, we obtain the additive renormahzation constant C (14.381) . 



^/ N (1 N + 2 Ur 



n — 1 ^ Cr 



127r 87r2 



(4.71) 



which gives the 0{u) correction to our previous result (I4.57p . Substituting into eq. (14.401) . 
we obtain 



Comparing the above result to eq. (14.581) . we observe that B receives no additional contri- 
butions at 0{u). Thus, from eq. (I4.62p . 



oiD NjN + 8) ^ c*. 



(4.73) 



6e V 127r 87r2 

which, upon determination of c* to order e would yield the entanglement entropy, eq. (14.461) . 



4- P{cr) to order 







a) 



b) 



c) 



d) 



e) 



FIG. 8: Contributions to {(fP'{x))n — {4''^{x))i at order (diagrams involving insertions of are 
not shown). The counterterm for the coupling u is shown as a circled dot. 



To complete our calculation, we need the value of the fixed point coupling c* to order 
e. This requires the knowledge of (3{cr) to order v?. As before, we will determine the 
renormahzation of c by computing the expectation value ([0^]r)n — ([0^]r)i- As explained in 
section IIV C 2\ we need to find only the additive renormahzation of c. Hence, we ignore all 
diagrams with vertices proportional to c^. At order m^, we obtain the graphs shown in Fig. 

El 
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Now we are faced with a new technical difficulty. Up to this point, to linear order in 
n — 1, the conical singularity entered our calculations through the tadpole term Gn{x, x) — 
Gi{x, x), whose form was fixed by dimensional analysis, eq. fl4.64p . up to an overall constant 
J{D). Moreover, the renormalization constants only depended on J{D = 4), which could be 
extracted from the explicit form of the propagator, eq. ( 14. 3p . However, at the present order, 
we are faced with the diagram in Fig. [8] a), which requires the full position dependence of 
the propagator Gn{x, x'). Yet, as far as we know, there is no simple expression for G„(a;, x') 
in arbitrary dimension, and even in D = 4 eq. (14.31) is rather awkward to work with. 

To address this problem, we expand the propagator Gn{x^x') to linear order in n — 1 
in terms of the usual propagators x'), eq. (14.661) . The simplest way to do this is to 

consider the 0{N) model in the presence of an arbitrary metric g^j,^, 

S = J rf^x^d^ + + ^0^^ (4.74) 

It is convenient to parameterize the n-sheeted Riemann surface using rescaled variables, 

f = ^/nr, if = 9/n (4.75) 

Then, the angular variable ~ + 27r. We may also define, 

f = r cos (fi, X = f sin (p (4.76) 

The coordinates (f , x) form the usual two dimensional Euclidean plane and uniquely specify 
each point on the Riemann surface. With this choice of variables, the metric (12.41) in the x\\ 
plane becomes, 

9a(3 = nSafs + - ^) (4-77) 
where a,f3 run over f, x. Note that we have chosen to rescale r in such a way that, 

det^ = l (4.78) 

Moreover, expanding g in powers of n — 1, g^/s = Sa/s + Sga/3, 

Sgap ^{n-1) (^6^^ - ^) (4.79) 

We drop the tildes on variables r, x in what follows. We can now obtain the usual Feynman 
graph expansion for the theory (I4.74p . treating 6gai3 as a perturbation. Note that all the 
integrals in the resulting expansion are over the usual D-dimensional Euclidean space. In 
particular, note that the bare propagator becomes, 

Gn{x, x') ^ Gi{x, x') + 6Gn{x, x') (4.80) 



24 



5G^{x,x') = {n-l) j d^'y \^6^^ - ^y^j d^G^ix - y)dpGr{x' - y) (4.81) 

By performing the integral, we immediately obtain eq. fl4.65p for x) — Gi{x,x). 

Using the expansion ( 14.8ip we compute the divergent part of the diagrams in Fig. [8] 
to linear order in n — 1. After accounting for the multiplicative renormalization of the 0^ 
operator, eq. (14. 48 p . we extract the additive renormalization of the coupling constant c, eq. 
fOTj) to 0(^2), 

^_n-l/ (iV + 2K , (iV + 2)(iV + 5) 7(iV + 2) n,^ ^^ 

12n [ e + e2 Sn^ 4e J ^^"^^^ 

and from eq. (I4.52p . 

P(Cr) = {N + 2)Ur ( 1 - ^ + (iV + 2)^ f 1 - Cr (4.83) 

Hence, 

27r ul\ , 27r 



_ 1 _ (,^ _ 1) = 1 _ (n - 1) (4.84) 

and from eq. (14. 73 p . 

N 

As = - — [n-l) (4.85) 

which by eq. (I4.46P finally yields the coefficient of the correlation length correction to the 
entanglement entropy, 

^ (4.86) 



1447r 



V. 4 - e EXPANSION: FINITE SIZE CORRECTION 



In this section we compute the geometric corrections 7, 7„ to the entanglement entropy 
and the Renyi entropy, eqs. (11.61) . (I1.14p . at the critical point. 

As before, we consider two semi-infinite regions A and B with a boundary at x = 0. 
However, we now take the remaining D — 2 spatial directions to have a finite length L. In 
order to avoid dealing with the zero mode, we use twisted boundary conditions along these 
directions. 

0(x + Mi) = e^'^»0(x) (5.1) 

where fii are unit vectors along the boundary. If the fields are real, then yjj = or tt. On 
the other hand, in an 0{N) model with N even, we can group our fields into A^/2 complex 
pairs - then, an arbitrary twist is allowed (however, this breaks the 0{N) symmetry down 
to U{1) X SU {N/2)). We note that when accessing D = 3 via 4 — e expansion, we will choose 
all ipi's to be equal. 

Thus, the boundary between regions A and B is a. D — 2 dimensional torus. Since this 
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manifold is smooth we expect the constants 7, 7„ to be universal. Moreover, we don't have 
to take into account divergences which appear as -D ^ 4 when the boundary has a finite 
curvature/ since this manifold is fiat. 



A. Gaussian theory 

Let us begin with the free theory. We wish to compute, 
7 N 

log^ = --iTr\og{-d\-nTr\ogi-d'),) (5.2) 

= - f E [^ni log(-5|| + kl)n - riT^ \og{-dl + kl),] (5.3) 

where k\ = ^2111^ and rii are integers. We leave the regularization of eq. (15. 3p implicit for 
now (we will later use dimensional regularization). Eq. (15. Sp involves the partition function 
of the two-dimensional massive gaussian theory evaluated in Ref. 0, 



log — 



= — (Trii log(— + m^) — TT-Trn log(— c?|? + m^)i) = — (n ) log(m^) (5.4) 

D=2 2 " " 24 V n 



Thus, 

'»S § = ^ (" - ^) E = (« - I) (X, X) (5.5) 



Here, Gj^{x,x) is the free propagator on an n-sheeted Riemann surface, which incorporates 
the finite size effects in the transverse direction. Explicitly, 

G^ni^,x') = J]G^=2(a:|i,4;fcl)e^^^-(-^--^^x) (5.6) 

In particular, for n = 1, 

justifying the last step in eq. (15.51) . 

An alternative representation for the propagator ( 15. 6p on the torus can be obtained by 
Poisson resumming k±, which is equivalent to "periodizing" the infinite volume propagator, 

G^{x, x') = J2 e^^^Gnix + IL, x') (5.8) 

I 

where / is a vector oi D — 2 integers in the plane parallel to the boundary. Note that when 
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X = x', only the / = term in eq. (15.80 is ultra-violet divergent and G^{x,x) — Gn{x,x) is 
finite. Moreover, since the / = term, Gi{x,x) ~ A^~^, is L independent, it gives a non- 
universal contribution to log(Z„/Z"), eq. (15. 5p . proportional to the area of the boundary. 
Concentrating on the universal constant term. 



log|^ = -N^ - - G^{x,x)) 

where from eqs. fl4.66p . (15. 8p . 



(5^9) 



L--(Gf(0)-G.(0)) = £<^^5:4;;i (5,10) 

Here and below we abbreviate x) by G^{0). 

We can now explicitly evaluate the universal constant contribution 7„ to the entanglement 
entropy for D = 3 and D = 4. 

In = fl + -] log(2| sm^/2\), D = 3 (5.11) 



12 V 

N 

7 = --log(2|sin(^/2|), D = ?> (5.12) 

D 

For D = 4, we note that the sum 

Y.C = {2nfG''=\^) (5.13) 

where G^^'^{0) is the massless two-dimensional propagator (with the zero-mode removed) on 
a torus with side-length 27r. This propagator can be expressed in terms of the Jacobi-theta 
function 9i, 

G--\0) = (log (^^,0 I - g -log^(o) (5.14) 

where rj is the Dedekind-eta function. 
Thus, 

ln = '^U + -\G^=\0l D = A (5.15) 



6 \ n 

^ = !?^G^=2(^)^ D = A (5.16) 



27 



B. 4 — e expansion 



We now compute the universal finite size correction to leading order in 4 — e expansion. 
The leading correction to the free theory behaviour comes from the boundary perturbation 
( l4.1Up . as at the fixed point c* ~ ^/e for 1 — n 3> e and c* ~ (n — 1) for |1 — n| -C e. Thus, 



^'''^og^ = -i] d^-'W{<P\r = 0)U = -^L--'Gt{r = r' = 0) = --^L^-^G{{x,x) 

(5.17) 

where in the last step we've used eqs. f l4.12p . (15.81) . Again, subtracting the non- universal 
area law piece ~ L^^^Gi(O), and combining eq. (I5.17p with the free theory result (15. 9p . 



Now replacing by it's fixed point value and taking D — > 4, 

Here we've set all the twists (pi equal. For 1 — n ^ e, eq. (14.240 . the c* term gives a correction 
of order ^/e to the free theory result. However, in the limit |1 — ^ e, eq. (I4.25p . the 
correction due to the boundary perturbation cancels with the free theory result to leading 
order in — 1, leaving. 

This implies that at the Wilson-Fisher fixed point the universal finite size correction to the 
entanglement entropy, 

7 ~ 0(e) (5.21) 
parametrically smaller than at the Gaussian fixed point in D = 4 — e. 



C. Beyond the leading order in e 

We now evaluate the universal finite size correction to the entanglement entropy 7 to 
order e. As before, we only work to leading order in n — 1. To order u the partition function 
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receives contributions from the diagrams in Fig. [91 The diagram in Fig. [9] a) is given by, 



c-0,1 1 

' log — 



4 



[{G^{x,x) - Gi{x,x)) + {G^{x,x) - Gi{x,x))] 
"~ - ^^^+/^"'-^\ Gf(0)-Gi(0))X: / d'x {G^=\x,x-kl)-G?=\x,x;kl)) 



N{N + 2){n - l)urfi' 
12 



{G^{o)-G^mj2w 

-. K I 



(5.22) 



where in the last step we've used eq. f l4.54p . 

The diagram in Fig. [9]b) can be evaluated with n = 1 propagators, 



6 ' log — 



(G't(O)- 6-1(0)) / d^-V / rfVGf(xx,x') 



'\2 



(5.23) 



Finally, the diagram in Fig. [9] c) can be obtained from eq. fl5.17p by substituting the 
counterterm for c, eq. fl4.68p . Combining all the diagrams in Fig. [9] with the 0(1) result, 
eq. (1231, 



log 



N /2 



(n-l) + c, )L^-2(Gf(0)-Gi(0)) 



(A^ + 2)Ur 

An 



{Lk^Y 2-ne 



(5.24) 
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Applying the usual technique for analytically continuing sums over D-dimensional vectors, 



1 r 

fc_L 

where 

n 

The function T{s) has the following asymptotics, 

T{s) s ^ (5.27) 

T{s) e-''p\ s oo (5.28) 

Hence, for finite (/? the integral in eq. fl5.25p converges in the s — > cxd region. Moreover, the 
s — >■ region contributes a pole for 4, 

E- — ; — -r h finite terms (5.29) 

As expected, this pole precisely cancels with the c counterterms, so that the expression 
f l5.24p is finite. Moreover, setting Cr to its fixed point value, eq. f l4.84p . the prefactor in eq. 
(I5.24P is already 0(e), so that we can neglect the 0{u) terms in the square brackets. Thus, 

and 

Note that the result fl5.3ip is of 0(1) in for A^ oo, instead of the naively expected 
0{N). It is not clear if this is an artifact of working to leading order in e. 

The function G^^'^{ip,(p) which determines the ip dependence of 7 is shown in Fig. [TOl 
We observe that 7 is a monotonically decreasing function of (p for < ip < it. In particular, 
for if = 71, 

Thus, 7 is negative for anti-periodic boundary conditions. On the other hand, for 9? -h> 0, 

A^e 

7~ r log 02, (j9 (5.33) 

' 6(A^ + 8) ^ ^ ^ 

suggesting that 7 is positive for periodic boundary conditions. Note that our expression 
for 7 becomes invalid for sufficiently small. The value of where the breakdown of 
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FIG. 10: The function ^((/j, if) determining the dependence of 7 on the twist 93, eq. (j5.3ip . 



direct perturbative expansion occurs can be estimated as foUows. Let us separate out the 
quasi- zero mode 0o of the field 0, 

'^(^) = JJ^U^y^''^"' + ^i^) (5-34) 

where (f){x) has the fc± = f mode omitted. At the mean-field level, the effective action for 0o 
is a two dimensional 0^ field theory, with an effective mass mlj-) ~ ^ and quartic coupling 
U2D ~ iD-i ■ We know that perturbative expansion in a 2D theory is valid for U2d/^2D 1- 
Thus, setting D = 4 and u = u*, we obtain, 

^p^ > e (5.35) 

as the domain of validity of perturbation theory. For smaller values of ip, the zero mode 
must be treated separately and non-perturbatively. This result can be checked in the 1/N 
expansion where one obtains a slightly stronger condition y?^ 3> e log 99. Cutting off the 
logarithmic divergence of (15.331) at the value of if where perturbation theory breaks down, 
we obtain. 

We conjecture that eq. (I5.36P is the leading order result for the case of zero twist (periodic 
boundary conditions). 
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VI. LARGE N LIMIT 



In this section we compute the correlation length correction to the Renyi entropy S'„, eq. 
f ll.lSp . in the large N limit. Although we are mainly interested in the physical case D = 3, 
we will keep the dimension of space-time arbitrary in our discussion in order to compare the 
results of the large- and 4 — e expansions. 

When working in the large- limit, it is more convenient to use the non-linear a- model 
version of the 0{N) model (12. ip . where the quartic interaction is replaced by a local con- 
straint (f)'^{x) = K Enforcing this constraint with the help of the Lagrange multiplier A(x), 
the action takes the form, 

S = J d^x(^^{d,<f>r + Ux{<j>'--^)^ (6.1) 

Our discussion in section IIIII is then directly transcribed into the present case with the 
replacement, t — (^ — y), 0^ iX. In particular, to determine the coefficient r„ of 
the correlation length correction to leading order in we need to find the behaviour of 
{iX(x)) at the critical point. 

We tune the 0{N) model to criticality g = Qc- Kt N = oo, the problem is reduced to 
finding the saddle-point value of the Lagrange multiplier {iX{x))n such that the gap equation, 

G.(x,x) = ^(02(a:))„ = -^ (6.2) 

is satisfied. Here G„(a;, x') is the Green's function of the operator —d"^ + {iX{x))n on the 
n-sheeted Riemann surface. The quantity Gn{x,x) requires regularization; we will implicitly 
use point splitting regularization. It is convenient to rewrite the gap equation as, 

Gn{x, x) — Gi{x, x) = (6.3) 

We note that at = cxd the scaling dimension of X{x) is 2, so, 

{zX{x))n = ^ (6.4) 

From (13. 4p . with the appropriate replacement 0^ —>■ iX, t —>■ g^^ — g^^, the constant a„ is 
related to the constant dn (13. 5p as 

dn = n_o ( — - - ] an (6.5) 



\9c 9 



Now from the gap equation at finite m, 

d^p f 1 1\ 1 



A^^ Ng^ J (27r)^ + J {An)^/ 



;T{1- D/2)m'^-'^ (6.6) 
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and 

,T{l-D/2)an (6.7) 



(47r)^/2 



In particular in D = 3, (i„ = -^a-n- Thus, the problem of computing the entanglement 
entropy &i N = oo reduces to finding the constants a„. 

We now need to find the Green's function G„. The main observation is that the angular 
harmonics on an n-sheeted Riemann surface are ^ — ^tie/n ^j^g^-g / jg integer. Hence, 

V27rn 

x') = I ^^e^^-(--^x)G^=2(r, r', O; kl) (6.8) 

where the two-dimensional massive propagator on an n-sheeted Riemann surface is given 
by, 

ii{e-e')/n 

G^=\r, r', e- m') = % 9i{r, r'; m') (6.9) 

I 

Here, 



r- 



We use spectral decomposition for gi, 



+ m \ gi{r,r';m ) = -6{r — r') (6.10) 



gi{r,r'-m^) = / f/E^^0,^^(r)0*^(r') (6.11) 



where 

ld_ fd_\{l/nl + a^ 
r dr \ dr j r 



riz ) + V 1 '^'-^ = E<P,E (6.12) 



and (pi^E are normalized to 



j drrct>lE{r)(t>i,E'{r) = 6{E - E') (6.13) 



The constant a„ must be positive in order to avoid the presence of negative energy states, 
which would render our saddle point unstable. Let us call the quantity I'^/n^ + a„ = u"^. Eq. 
f l6.12p admits two linearly independent solutions. 



J\u\iVEr) (6.14) 



0(r) = i=J_|H(v^r) (6.15) 

We recall that 

J^(x) ~ |xr, X (6.16) 

When working in free space (in the absence of conical singularity and potential (16. 4p ) one 
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chooses only the solutions with a positive index |z/| = |/|, so that is finite and differ- 

entiable at r = 0. However, in the present problem there is no a priori physical reason why 
the solutions (and hence the propagator) have to remain finite as r ^ 0. 

In fact, a particle in a potential is a famous problem known as conformal quantum 
mechanics. Note that the potential ( 16. 4p is highly singular and requires regularization at 
short distances. Such regularization will automatically appear in the linear 0{N) model, 
which can be obtained from (16. ip by adding a term A^/4u to the Lagrangian. In that case, 
eq. (16.41) only holds for ur^~^ ^ 1 and the saddle point value {i\) is modified at short 
distances. We note that even after this regularization, the / ^ states still experience an 
/^/r^ centrifugal barrier and we must choose positive index solutions (16.141) for them. We 
now concentrate on the / = sector. For simplicity, imagine cutting the 1/r^ divergence 
off at some radius r = tq and replacing it by a finite potential. Generally, the resulting 
scattering states will approach the positive index solutions (I6.14p for VEtq 0. However, 
non-trivial behaviour can occur if the potential is close to developing a bound state. In that 
case, for < 1, one "dynamically" generates a length-scale ^ and the scattering solutions 
become linear combinations of (I6.14p and (16.150 with coefficients (and, thus, the phase-shifts) 
depending on VE^. Since we are looking for a scale invariant solution to the gap equation, 
we need ^ oo, i.e. the system is exactly at the threshold of bound state formation. At 
this threshold, for \/ErQ one obtains negative index solutions (I6.15p . Note, that this 
behaviour is special to the range < 1 and does not occur for |z/| > 1. This fact could 
be anticipated as the negative index solutions are square integrable at short distances for 
|z/| < 1 but not for \i>\ > 1. 

Thus, applying RG terminology to the simple quantum mechanics problem ( I6.12p . we 
conclude that there are two fixed points - one stable (I6.14p and one unstable (I6.15p . However, 
we are allowed to choose the unstable fixed point solutions as we are fine tuning both the 
long and short distance parts of {iX) to solve the gap equation. 

With these remarks in mind, 

/■°° 1 

gi{r,r';m^) = / kdk— -J,^{kr)J,Xkr') (6.17) 

where vi = a for / = and vi = ^yp/n"^ + for |/| > 0, with = a^. The constant a can 



be either positive or negative. We note that as discussed in Ref |l3|, a enters the operator 



product expansion of the field 0( approaches the conical singularity, 

<j){xi\,x±) ~ rXO, x±), r (6.18) 
Combining eqs. (16. 81) . (16. 9p and (16.171) . and performing the integrals over k±, k we obtain 
G ir-r'ex, -x')- ^((3-/^)72) r(D/2 - 1 + 

Since we are mostly interested in Gn{x = x'), we have set r = r' , x± = x'j_ in (I6.19P : we have 
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left ^ 7^ as a regulator. 

As an aside that will be of some interest later, we note that fl6.19p is meaningful only for 
a > —{D /2 — 1). For a<— (Z)/2 — 1) one obtains an infrared divergence in the /c^, — > 
region of integrals (16.81) . ( pTTTj) . We note that at a = —{D/2 — 1), eq. (16.121) has a zero 
energy solution, 

= ^ (6.20) 

The solution (I6.20p could, in principle, correspond to a saddle point with a non-zero ex- 
pectation value {(f){x)). Note that the r dependence of (I6.20p is consistent with the scaling 
dimension [0(x)] = D/2 — 1 in the ^ oo limit. Alternatively, observe that the scaling 
dimension of the "boundary" operator, [0(0, x_l)] = D/2 — l-|-a— >-Oasa— > —{D/2 — 1), 
indicating a tendency to condense. However, the infrared divergence of the propagator 
(I6.19P indicates that condensation of 0(x) at the conical singularity is unstable to fluctua- 
tions. This is not unexpected, as the condensate would he D — 2 < 2 dimensional. Such 
a condensate certainly cannot exist for any g > Qc, a.s it would violate the Mermin- Wagner 
theorem. Long range interactions could potentially stabilize the condensate exactly at the 
critical point, however, the above discussion shows that this does not occur (at least in the 
large- limit). 

We use contour integration to write (16.190 in a somewhat more convenient form, 

Gn{r = r',e,x± = x'_l) 

(6.21) 



with 



T{-iu + D/2-1) T{iu + D/2-1] 



cosh(i/(7™ — \6\)) 
smh.{7mu) 

-iT{3 - D) ^ 2 - D/2) T{iv + 2 - D/2)_ 

27rr(3-D)sin(7r(3-L')/2)sinh(7rz/) 1 

coshVz/-sin2(7r(3-L')/2) \V{iu + 2 - D /2)\'^ 



(6.22) 
(6.23) 
(6.24) 



In particular, for = 3, -R(z^) = 7rtanh(7rz/). We note that despite the presence of the 
0{—a) term in eq. (16.2 ip . Gnir = r' , 9, x± = x'j_) is analytic at a = as is evident from eq. 
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fl6.19p . Thus, the gap equation (16. 3p takes the form, 



(Gn — G 



I x = x' 



47r^/2r(2-D/2)r^-2 



+ 9(-a)-R{ia 
n 



du ( = coth(7mV z/2 + a^) — coth(7rz/) ] R{u) 



(6.25) 



The function -R(z/) is positive for real values of v. So the left-handside of the gap equation 
goes to —OO as a — > OO and to cxd as a — > —{D/2 — l)"*". Hence, the gap equation always has 
at least one solution, and more generally, an odd number of solutions. Numerically, we find 
that the gap equation has a unique solution for all n ioi D < D^, ~ 3.74. For D > D^, 
there are one or three solutions depending on the value of n, as we will discuss below. 

As we are mainly interested in the entanglement entropy, let us consider the limit n 1. 
Then we expect a 0. The integral in (16.251) is non-analytic at a = 0, due to singular 
behaviour in the z/ — > region. Noting that -R(z^) ~ -R'(0)z/, as z/ — > 0, we obtain to leading 
order in a, 



{Gn — Gi)\ , — {Gn — G 



x=x' ,a=0 

- Q{—OL)a 



R'{o) (I r, f y'' . 



4n7r^/2r(2 - D/2)r^-2 V 2' ' ' 'J An^/^r{D - l)r^-^ n 

(6.26) 

where the contributions from the integral and the 6 function have combined to produce a 
result analytic in a. Now using eqs. (14.641) . (I4.65P for {Gn — Gi] 



I x=x'a=0 ' 



n^l (6.27) 

Note that the exponent a controlling the OPE (16.180 of the field (^{x) at the conical singu- 
larity is positive for n < 1 and negative for n > 1. Now, from (16.270 . 



{D-2) 



2 



(n-l)2, n^l (6.28) 



" A{D - 1)2 
Therefore, combining eqs. (13. lip and (16. 5p . we find that 

Tn OC — OC 77. — 1, n ^ 1 (6.29) 

1 — n 

and the correlation length correction to the entanglement entropy proper vanishes at leading 
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order in 



r = lim r„ = (6.30) 

Thus, for all dimensions 2 < D < A 

r ~ 0{N) (6.31) 

even though r„ ~ 0{N'^) for all n 7^ 1. 

So far we have concentrated on the solution to the gap equation in the n 1 limit for 
arbitrary dimension. However, we can also obtain an analytic solution for arbitrary n in the 
limit D = 4 — e. Such a solution is useful for comparison to the results of the 4 — e expansion 
presented in section HVl 

When D = A — e, the function -R(z/) = — 2z/^~'^r(— 1 + e). The divergence of the F function 
is not important here as it is just an overall factor in the gap equation (which anyway 
cancels with r(2 — D/2) in (16.211) ). However, the integral (16.251) now diverges for z/ — 00 
if e = 0. Hence, for generic n and D = 4 — e the leading a-dependent contribution to the 
gap equation comes from the region z/ ^ 1 and is of order, -^a^. This suggest that a will 
be at most of order e^^^. However, for a very small (i.e. n — 1), we already know from 
the previous discussion that the leading contribution to the integral scales as |a| and comes 
from the z/ — > region. Keeping these two contributions (one non-analytic in a and the 
other analytic, but with a diverging coefficient) and setting a = in the rest of the integral, 
we reduce the gap equation to 



Jo \z/^ + a^ / Jq \sinh(7™z/) sinh(7rz/) 

1 /'OO 

- -a^ / du ly-'-^ - e(-a)- = (6.32) 

n 



n e 6 \n'^ 

The quadratic has two solutions 



a 1 / 1 , , , 

r + --J ITT-l =0 (6.33) 



and the corresponding values of dn, eq. (16. 7p . are. 



" 87r2 



6 Kn"^ J 2n^ 



(6.35) 



Eq. (16.351) is in agreement with the result of the 4 — e expansion, eq. (I4.27p . and we can 
identify the a± saddle points with the cf fixed points. Moreover, we see that the predictions 
of the large- (I6.34p and 4 — e expansion (14.340 for the OPE exponent a also agree. Note 
that both saddle points (16.341) disappear for n > ric ^ 1 + 3e/4. This coincides with the 
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n 


a„ 


2 


-0.16515 


3 


-0.26594 


4 


-0.32905 


5 


-0.36743 



TABLE I: Solution to the gap equation in the large- limit for D = 3. 



value of n at which runaway of RG flow is observed in the 4 — e expansion. However, as 
we noted earlier, the gap equation always has an odd number of solutions. Thus, we have 
missed a solution in our discussion above. This solution has a ~ —{D/2 — 1) ^ —1, i.e. a 
is not small. Its existence is possible due to a cancellation of 1/e divergences between the 
large u part of the integral and the d{—a) term in (16.251) . Keeping these two contributions 
to the gap equation, we obtain in the a —{D/2 — 1) limit. 



e na + D/2-1 

So, 



(6.36) 



a = -l + le + -e^ (6.37) 
2 n 

Eqs. (I6.34p . (I6.37P comprise the three solutions to the gap equation for 1 < n < ric, and 
eq. ( I6.37P is the only solution for n > ric- We speculate that the runaway of the RG flow 
observed in 4 — e expansion for n > ric is towards the flxed point (16.371) . As we noted 
above, the value a = —{D/2 — 1) corresponds to the would be condensation of the fleld 
at the conical singularity. Thus, for e — *• 0, the saddle-point (16.371) is proximate to such 
condensation. This is consistent with our interpretation of the RG flow c — — oo as the 
tendency to formation of {4>{x)) ^ 0. However, the large- analysis demonstrates that no 
true spontaneous symmetry breaking at the conical singularity occurs for D < 4. 

To our knowledge no such non-trivial n-dependence has been previously observed in any 
theories. Still, in the large- A^ expansion such behaviour is only present for > ^ 3.74 
and its relevance to the physical case D = 3 is doubtful. Moreover, the non-analyticity 
occurs away from n = \ and, thus, is unimportant for computing the entanglement entropy 
proper. Indeed, the behaviour of the theory for n — >■ 1 (I6.27P is found to evolve smoothly as 
the dimension D increases from 2 to 4. 

We now come back to the physical case D = 3, where the solution to the gap equation is 
unique. The numerical solution for the flrst few integers n is listed in Table [B Then, from 
fl3:T0|) and fl63D . 

37r2Ar2 nal 

The coefficient (16.381) can be, in principle, obtained numerically by performing classical 
Monte-Carlo simulations of the 0{N) model in the spirit of Ref. QJ. 
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So far our large- computation has been confined to tlie correlation length correction to 
the Renyi entropy. At leading order the calculation was technically fairly simple, as utilizing 
the discussion in section UTTI we could work at the critical point. In particular, the form of 
the Lagrange multiplier (zA(r)) was fixed by scale invariance up to an overall constant. To 
proceed beyond the leading order, as is required for the calculation of the correlation length 
correction to the entanglement entropy proper, we would have to work in the gapped phase. 
The Lagrange multiplier (zA(r)) would now be a non-trivial function of r with a length scale 
determined by the correlation length ^ = m~^. Similarly, if we wish to compute the finite 
size correction 7 to the entanglement entropy, (iA(r)) will again vary non-trivially with a 
length scale determined by the size L of the compact direction. In both cases, we have 
to solve the gap equation for a whole function (iA(r)) rather than a single number a„. In 
principle, this problem can be addressed numerically. It would be particularly interesting 
to check whether 7 ~ 0(1) for — >• 00 as suggested by the 4 — e expansion, eq. fl5.3ip . 



VII. CONCLUSION. FUTURE DIRECTIONS. 

In the present work we have computed the universal finite size and correlation length 
corrections to the entanglement entropy and the Renyi entropy for the 0{N) model. The 
evaluation of this entropy required a study of the 0{N) field theory on a n-sheeted Riemann 
surface for general n, and an understanding of the nature of the n ^ 1 limit. For n 7^ 1, 
there is a conical singularity at the origin of the Riemann surface and we have presented a 
detailed analysis of the structure of the "boundary" excitations of the 0{N) CFT at this 
singularity. (A closely related CFT with vortex boundary conditions was studied in Ref. Il3 
with a very different physical motivation.) In particular, we showed that in the context 
of e = 4 — D expansion, the RG flow of the boundary couphng c in Eq. (I4.10p was the 
key to a determination of the entanglement entropy. The RG flow of c had two possible 
structures shown in Figs. [5] c) and d). For n greater than a critical n^, we had flow in 
the infrared to c = —00 as in Fig. [5] d). In contrast for n < Uc, we had three possible 
fixed points, and the n ^ 1 limit was controlled by the non-zero fixed point c = c^, at 
which all strong hyperscaling assumptions were obeyed. All our computations in the e and 
expansions were consistent with this RG flow and flxed-point structure. One crucial 
consequence of the boundary perturbation and the subtle limit n — 1 is that the flnite size 
and correlation length corrections to the entanglement entropy are different at the Wilson- 
Fisher and Gaussian fixed points already at leading order in e expansion. 

In this paper we have considered a geometry with a smooth, straight boundary between 
regions A and B. One possible extension of our work is to consider boundaries with sharp 
corners. In such geometries, it is expected that the entanglement entropy will contain a 
universal logarithmically divergent term.-i^i^ Moreover, we have only studied the correlation 
length correction to the entanglement entropy in the symmetry unbroken region t > 0. It 
would be interesting to extend our treatment to the symmetry broken phase t < 0. 

While our paper was being completed, we learned of the numerical study of entanglement 
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entropy in the d = 2 quantum Ising model in Ref. Il5|. At the quantum critical point the 
authors of Ref. 
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find evidence for a finite size correction 7 as in Eq. (11.61) in the case when 
the boundary between regions A and B is smooth. We note that the geometry studied in 
Ref. Il5| is an L X L torus divided into two equal cylinders rather than the infinite cylinder 
cut in half that we have considered here. Thus, the two results cannot be compared directly. 
Nevertheless, the value of 7 in Ref. [l5| is found to be positive, as in our conjecture in Eq. (II. 8p 
for the case of periodic boundary conditions along the cylinder. 
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